Quantum-wave evolution in a step potential barrier 



(N 
O 

o 

» 

o 
O 



> 
O 

o 
o 

O . 
^ ' 

^— > ■ 
C ■ 

=5 

cr 



X 



Jorge Villavicencio, 1 ' 2 '0 Roberto Romo, 2 '[] and Sukey Sosa y Silva 2, [] 

1 Instituto de Fisica, Universidad National Autonoma de Mexico 
Apartado Postal 20 364, 01000 Mexico, Distrito Federal, Mexico 
2 Facultad de Ciencias, Universidad Autonoma de Baja California 
Apartado Postal 1880, 22800 Ensenada, Baja California, Mexico 
(Dated: February 1, 2008) 

By using an exact solution to the time-dependent Schrodinger equation with a point source initial 
condition, we investigate both the time and spatial dependence of quantum waves in a step potential 
barrier. We find that for a source with energy below the barrier height, and for distances larger than 
the penetration length, the probability density exhibits a forerunner associated with a non-tunneling 
process, which propagates in space at exactly the semiclassical group velocity. We show that the 
time of arrival of the maximum of the forerunner at a given fixed position inside the potential is 
exactly the traversal time, r. We also show that the spatial evolution of this transient pulse exhibits 
an invariant behavior under a rescaling process. This analytic property is used to characterize the 



evolution of the forerunner, and to analyze the role played by the time of arrival, 3 
recently by Muga and Biittiker [Phys. Rev. A 62, 023808 (2000)]. 
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I. INTRODUCTION 



Since the original proposal by Stevens |l| of tunneling 
monochromatic fronts, the problem has become contro- 
versial and still open to investigations. The existence of 
such propagating fronts has been supported in asymp- 
totic analysis based on exact analytical solutions to 
the time-dependent Schrodinger equation, with point 
source initial conditions. However, analytical and numer- 
ical studies based on different boundary conditions and 
potentials || , have shown evidence that it is not possible 
to identify such fronts, in the way suggested in Ref. 0. 
Despite these clarifying works, the need of a direct and 
comprehensive study of the original point source problem 
has been recently emphasized [Q, |E| . 

A recent study of Muga and Biittiker Q of quantum 
waves in a potential step, has shown the existence of 
a traveling transient forerunner for frequencies of the 
source below the barrier height, which in opaque bar- 
rier conditions (xqo 3> 1) is dominated by over-the- 
barrier frequencies, implying that this transient struc- 
ture is associated to a non-tunneling process. They also 
found that the time of arrival at a fixed position xj is 
unexpectedly given by 3 _1 / 2 r, where r is the traversal 
time, defined as r — (xf/v qa ), with v qo = (hqo/m), and 
go = [2m(Vo — Eq)] 1 / 2 /h; the parameters Vo and Eq cor- 
respond respectively to the step potential height and the 
initial energy of the source. We believe that in order to 
gain more insight about the evolution of quantum waves, 
the above investigation should be extended to study the 
role played by r in the spatial-dependence of the solution. 
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In this paper we investigate both the time and position 
dependence of quantum waves in a step potential bar- 
rier by means of an exact solution to the time-dependent 
Schrodinger equation. Our analysis shows that the fore- 
runner travels in space at exactly the semiclassical group 
velocity v qo , implying that the traversal time t is exactly 
the time of arrival of the maximum of the forerunner at 
a fixed point deep inside the potential step; clearly this 
result is different from the time scale of Ref. || by a fac- 
tor of 3 -1 / 2 . Using a scaling property of the forerunner 
derived from the analytic solution, the exact interpreta- 
tion of the time scales r and 3~ 1 ^ 2 r as times of arrival is 
clarified. It is a surprising fact that r, a time scale tradi- 
tionally regarded as one of the possible tunneling times, 
is here associated to a non-tunneling process. 

This work is organized as follows. Section II deals with 
numerical examples, showing both spatial and time evo- 
lution of the probability density for initial waves with 
energies Eq > Vq, and Eq < Vq. Finally, the conclusions 
are presented in Section III. 



II. QUANTUM WAVE EVOLUTION 

In this section we shall explore the main features of 
quantum wave evolution. Here we use the exact analyt- 
ical solutions ||] (see also Appendix A) to the time- 
dependent Schrodinger equation for a step potential bar- 
rier V(x) = Q(x)Vq, for a wave formed by a point source 
with a sharp onset, 



ipo (x = 0, t) 



-iujQt 

o, 



t > 0, 
t < 0. 



(1) 



We have defined loq = (Eq/H), where Eq corresponds 
to the initial energy of the source. It is worthwhile to 
point out that the source boundary condition given by 
Eq. (nj) is not a standard one in quantum mechanics. 
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For a physical interpretation of these type of boundary 
conditions see Ref. 

The analytical solutions for the point source problem 
are given by, 

ip>(x,t) =e- lVt [M( Xl k Ql t) + M{x,-k ,t)}, loo > V, 

(2) 

and, 

ip< (x, t) = e~ lVt [M(x, iq , t) + M (x, -iq , t)] , uj < V, 

(3) 

where fc = [2m(wo — V)/h] 1 / 2 , q — ik , and V = (Vo/h). 
In the above equations the M— functions are defined as, 



M(x, q, t) ee M(y q ) = l -e m ^l 2ht w{iy q ), (4) 



where w is the complex error function [ p"o[ |IT 
exp(— z 2 )erfc(— iz). The argument y q in Eq. 
by, 



Vq 



\2ht) 



hq 
x 1 



rn 



w(z) = 
) is given 



(5) 



where q — ±/co, ±z<Zo- 

From the analysis given in Subsection A.l, one can 
obtain the long time behavior of Eqs. (|^) and (||). For 
the long time regime t — > oo the asymptotic solutions 
are, 



and 



i>%(x,t)=e 



LO < V, 



LOq > V. 



(6) 



(7) 



In what follows we shall explore the main features of 
waves along the potential region, x > 0, generated by a 
source with the initial frequencies luq > V and too < V. 
Our emphasis will be in the second case, where we shall 
reveal novel aspects in the spatial evolution of quantum 
waves. 
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FIG. 1: (a) Spatial evolution of |t/j> | , for two different val- 
ues of time: ti — 15.0 fs, and t-z = 30.0 fs. The positions 
x sc (ti) = Vk ti of the semiclassical fronts are indicated by 
an arrow in each curve, (b) Time evolution of |'0>| 2 at 
different fixed positions x\ — 5.0 mn and X2 = 30.0 nm, 
in order to show the wavefront propagation. The positions 
t sc (xi) = (xi/vk ) of the semiclassical fronts are also indi- 
cated by an arrow. Notice how as time increases the solution 
tends to the correct asymptotic value given by \ip>\ 2 = 1. 



A. Source with an initial frequency luq > V 

Although the main purpose of the present work is to 
explore the dynamics of the solution for a source with an 
initial frequency loq < V, for completeness we begin our 
discussion by illustrating some cases for loq > V. In this 
subsection we investigate the evolution of the solution 
given by Eq. We choose the following parameters 

for the system: height of the potential step, Vq = 1.0 
eV, and the energy of the initial wave, Eq — 2.0 eV. The 
calculated values of |^>| 2 are shown in Fig. [I] (a) as a 
function of the position x, for two fixed times t\ = 15.0 
fs, and ti = 30.0 fs. It is clearly appreciated in Fig. 
(a), that a main wavefront with damped oscillations 
is propagating from left to right. One can easily verify 



by numerical inspection that the main wavefront trav- 
els approximately with the semiclassical group velocity 
Vk — hko/m. The arrows in Fig. [j] (a) indicate the 
semiclassical positions x sc {ti) — Vk ti,(i — 1,2). 

The plots of |V'>| 2 vs t, for a two different fixed po- 
sitions x\ — 5.0 nm and X2 = 30.0 exhibited in Fig. 
^ (b), complement our discussion for this case. Here 
we also illustrate the propagation of the main wave- 
front for the specific positions Xi indicated in the fig- 
ure. Notice how the probability density rises from zero 
and grows monotonically towards a maximum value, from 
which it starts to oscillate until it reaches the stationary 
value. This sudden rise occurs approximately at times 
t S c(%i) = i x i/ v k )i which correspond to the arrival of the 
main wavefront at the fixed position xi. It is interesting 
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to mention that these propagation features are similar 
to the diffraction in time patterns observed in free parti- 
cle propagation within the Moshinsky shutter model ||, 
which is in fact a different type of initial condition. 

Although we see that in the case u>q > V it is per- 
missible to speak of a propagation of a main wave front, 
the issue has become controversial in case of initial waves 
with uj < V 0, §, |, H @. The results for this more 
complicated case is discussed in the next subsection. 



B. Source with an initial frequency loq < V 

In what follows we shall explore the main features of 
the time and spatial behavior of the probability density 
for a source with an initial frequency ujq < V. We con- 
sider here the same potential parameters as in the previ- 
ous subsection, but the energy is now chosen below the 
height of the barrier, Eq — 0.5 eV. In Fig. || (a) we show 
the time evolution of |V'<| 2 at the positions xi — 1.2 nm 
and X2 — 1.5 nm. Notice that the curves exhibit a be- 
havior similar to the diffraction in time pattern observed 
in the previous subsection. As we shall discuss below, 
this behavior is not related to the propagation of a main 
wave front, and requires a different interpretation. In 
fact, the time-diffraction- like pattern disappears after a 
certain depth inside the potential, where the probabil- 
ity density exhibits now a pulse-like structure, as can be 
appreciated in Fig. || (b). 

In order to understand the transition between the quite 
different dynamical behavior of \ip<\ 2 observed in Fig. ^ 
(a) and fb), a series of snapshots of |V'<| 2 vs x are shown 
in Fig. @ (solid lines) for increasing values of time. We 
can appreciate in Fig. || (a) that at short times the prob- 
ability density behaves like the stationary solution |'i/'<| 2 
(dashed line) along the internal region. However, the 
probability density, |"0<| 2 j dramatically changes as time 
elapses. This is depicted in Fig. [| (b), where we see 
the appearance of a bump for values of the position x 
larger than Xq. This special point Xq corresponds to 
the onset of a pulse-like structure that we call transient 
pulse (forerunner). In Figs. |^ (c) and (d), the transient 
pulse is now fully formed, and propagates along the re- 
gion x > Xo. However, notice that as time increases, the 
birth of the pulse occurs at values of the position Xr, 
greater than value of the lower-bound, Xq. The posi- 
tions X (full square) and Xr (full dot) are indicated in 
Fig. ||. We have systematically studied this transition 
for different potential steps, as well as for different ener- 
gies of the initial waves, and our results indicate that the 
onset of the transient pulse always occurs at a specific 
depth, Xq, inside the potential, which is twice the value 
of the characteristic penetration length x p = Qq 1 of the 
stationary solution: Xq m 2x p . The existence of this two 
regimes for the point source model has also been reported 
in Ref. H. A similar separation into two regimes, but in 
different models has been observed by one of the authors 
inRefs. @,|l|. 
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FIG. 2: Time evolution of |V , <| 2 ) at different values of the po- 
sition along the internal region. Notice that (a) at the fixed 
positions x\ = 1.2 nm, and X2 = 1.5 nm, the probability den- 
sity exhibits an oscillating behavior similar to the diffraction 
in time phenomenon. At the positions (b) X3 = 6.0 nm, and 
24 = 10.0 nm, located deep inside the internal region, the 
oscillating pattern no longer exists, and it is replaced by a 
pulse-like structure. 

We shall later see that the approximate value of Xr 
can be determined approximately by using a simple ana- 
lytical formula. In what follows we shall explore in more 
detail the regimes < x < Xq, and x > Xq. We shall 
refer to the latter as the transient pulse regime. 

1. < x < Xo regime 

This regime corresponds to short distances inside the 
potential, namely, to positions x such that < x < Xq. 
Here we shall explain the oscillatory patterns of the \ip < | 2 
vs t plots, observed within this finite interval and ex- 
hibited in Fig. | (a). We shall see that, rather than 
originated by the passage of a traveling wavefront (as it 
occurs in the case ujq > V), this effect is due to fluc- 
tuations of the probability density. This dynamical be- 
havior is illustrated in Fig. |], where we present a se- 
ries of snapshots of |V'<| 2 vs x, for increasing values of 
time indicated in the figure. Here we have included for 
comparison the calculation of the stationary probability 
density |"0<| 2 (dashed line). Notice that the evolution of 
the time dependent probability density occurs in such a 
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FIG. 4: Plot of \ip< | as a function of the position at different 
values of time: to = 1.0 fs (solid line) ti = 2.0 fs (dotted 
line), and t% — 3.0 fs (dashed-dotted line). The stationary 
solution |V><| 2 (dashed line) is also included for comparison. 
The inset shows \ip<\ 2 (solid line) at a later time tz = 15.0 fs, 
which becomes indistinguishable from the stationary solution 
(dashed line). 



way that |f/ ; <| 2 is found sometimes above or below the 
asymptotic value. This illustrates that, in this regime, 
the probability density fluctuates around the stationary 
solution before reaching its asymptotic limit. At long 
enough times, we must have J'0<| 2 — ► |i/'<| 2 , which is il- 
lustrated in the inset of Fig. H, where for t = 15.0 fs both 
curves become indistinguishable. 

The above results are relevant since they clearly show 
that we can not associate in this regime a semiclassical 
group velocity to the main wave-front towards its steady 
state. In other words it is not possible to define a semi- 
classical main front attenuated exponentially by e~ q ° x , 
traveling along the internal region with a group velocity 
v qo = (Hqo/m), as suggested by Moretti ||. However, 
a special kind of propagation actually exists for x > Xq 
(opaque barrier conditions), which consists of the travel- 
ing forerunner mentioned in the Section I. It should be 
emphasized however that this transient structure (in view 
of its frequency content) is not associated to a tunneling 
process ||. We shall show in the next subsection that the 
maximum of this forerunner travels at exactly the group 
velocity, v qo . 



2. Transient pulse regime (x > Xq ) 



FIG. 3: Birth of the probability density |i/'<| 2 (solid line) 
as function of the distance x, for increasing values of time: 
ti = 1.0 fs, t a = 3.0 fs, tz = 4.0 fs, and U = 15.0 fs. In part 
(a), we see a monotonically decreasing curve. However, a 
bump begins to appear in (b), which corresponds to the birth 
of a transient pulse. In (c) and (d), the maximum of this 
pulse propagates along x > Xq. In all cases, the stationary 
solution \tp< 2 (dashed line) is included for comparison. A full 
square and a dot indicate the position of Xq = 2.134 and Xr, 
respectively. 



This regime corresponds to the case x > Xq, where 
the dynamics of the probability density \ip < | 2 is governed 
by a transient pulse. This transient structure originates 
from Xq onwards, and evolves thereafter along x > Xq, 
as depicted in figures | (c) and | (d). In Fig. | (a), we 
plot \tp < | 2 vs x for the increasing values of time: to = 30.0 
fs, h = 100.0 fs, t 2 = 150.0 fs, and t 3 = 300.0 fs. Here, no 
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FIG. 5: (a) Spatial evolution of the transient pulse calculated 
from \tp<\ 2 for increasing values of time: io = 30.0 fs, ti = 
100.0 fs, i a = 150.0 fs, and t 3 = 300.0 fs. Notice the different 
values of the pulse peak and width for each particular time, 
(b) Plot of the rescaled probability density rii\tp < \ 2 (i = 1, 2, 
3) as function of a variable x/rji, for the same cases depicted 
in (a) . All curves become indistinguishable and coincide with 
the \ip<\ 2 vs x plot of the case t = to- Also included here is 
a plot of Eq. ([L3j), which, as expected, coincides exactly with 
all the curves. 



|V><| 2 vs x corresponding to to- 

The scaling property illustrated above by means of a 
numerical example, can also be explained using our an- 
alytic formulas. However, instead of using the formal 
solution given by Eq. (^), we shall first derive a simple 
expression for the spatial evolution of |V'<| 2 - Notice that 
the behavior of |V><| 2 is essentially governed by two M 
functions, namely, M(yi qo ) and M(y_j 9o ). By a numeri- 
cal inspection we find that the series of the M functions 
discussed in Subsection 1 of Appendix A, apply in this 
regime. In order to exploit the analytical properties of 
the Moshinsky functions, first we have to determine the 
location in the complex-plane of the phases 4>±i qo of the 
arguments y±i qo ■ We find that the phase <\>-i qa , is always 
within the interval — n/2 < <fi-i qo < tt/ 2. Thus the se- 
ries representation given by Eq. (AlOa) applies in this 
case. Similarly, we find that although the phase 4>i qa , may 
change from tt/2 < (j) iqo < 3n/2 to — tt/2 < (j) iqo < n/2 
in the vicinity of the maximum of the transient pulse, 
t he ser ies representation of the M functions given by Eq. 
(Alia) gives an excellent approximation. Therefore, we 
find that in this regime, the behavior of the probability 
density is essentially governed by the first terms of the 
series representation of Eqs. (AlOa) and (Alia). That 



and 



M(y-i 



M(y lqo ) ~ -e™ 



2 /2ht 



'/2ht 



2e Viq o 



(8) 



(9) 



After substituting the above expressions in Eq. (g), and 
performing a few algebraical manipulations, we obtain 
the following expression for the probability density, 



4> < ~ v< + iPt P , 



(10) 



where ip< is the stationary solution given by Eq. (||), 
which decays exponentially along the internal region. 
Here iptp stands for, 



traveling main front is observed, as in the ojq > V case. 
Instead, we see a forerunner that broadens as it travels 
along the internal region, with diminishing amplitude, 
similar to the propagation and spreading of a wavepacket. 

In addition to the above, the evolution of the proba- 
bility density in this regime, has an interesting scaling 
property. If we express the times ti (i = 1, 2, 3) as multi- 
ples of £o j namely ti — rjito (in our examples: r]i — 3.333, 
r/2 = 5.0, and 773 = 10.0), and use the factors rji to renor- 
malize both axis, in such a way that we now plot 7?i|V'<| 2 
vs (x/rji), a striking result is obtained: all curves coin- 
cide. This invariance under scaling is illustrated in Fig. |5| 
(b), where all the curves corresponding to i — 1, 2, 3, are 
indistinguishable among them, coinciding with the curve 



1 



tp 



2tt 1 /2 



i(mx 2 /2ht-Vt) 



1 

Viqo 



1 



y-i 



(11) 



which as we shall see below, accurately describes the tran- 
sient pulse. We have found by numerical inspection that 
Eq. ( |l0| ) gives an excellent description of the spatial evo- 
lution of |^<| 2 , provided that t > (Xo/v qo ). The result 
given by Eq. (ITy) shows us that the probability density 
can be described by a subtle interplay between -0" and 
iptp- A measure of such an interplay can be given by 
the ratio R = /ifj t p\ 2 - We have found that the ob- 
served transition discussed in Fig. |], occurring at Xr, 
corresponds to values of the position such that condition 
R = 1 is satisfied. 

If we now focus our attention in describing the dynam- 
ics of the pulse, from Eq. (Ill]) we can obtain a simple 
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FIG. 6: Comparison of l^tp] 2 (solid line) and |V'<| 2 (dashed 
line), (a) as a function of time, for xq = 8.0 nm, and (b) as a 
function of position, for to — 30.0 fs; in this case, both curves 
overlap, and become indistinguishable among them. 



analytical expression for this transient structure, namely, 



|^ tp (x,t)| 2 = - 



(hx 2 t/m) 



1" [x 2 + (Hq t/m) 2 



(12) 



The analysis from Eq. (|3|) to Eq. (|12|), is equivalent to 
the one performed in Ref. || in terms of a pole and a 
saddle-point contribution, with a different notation and 
terminology. It can be appreciated in Fig. || (a), that 
Eq. ( [12] ) provides an excellent description for the propa- 
gating pulse. Here, we have plotted iV'tpl 2 (solid line) as 
a function of time at a fixed position x — 8.0 nm, and 
the curve almost coincide with the calculated with the 
formal solution, Eq. (^) (dashed line). We also compare 
Eqs. (|J) and ( |l2| ) with a plot of the probability density 
as a function of position, see Fig. |](b), and an the de- 
scription of our approximate formula is excellent, since 
both calculations are indistinguishable in the graph. 
In order to show the scaling property, we make in Eq. 



( |l2| ) the general substitution of t' — rjt and x' — r\x. The 
result of this rescaling is the following equation, 



(hx 2 t/m) 



W [x 2 + {hq t/m) 



(13) 



A plot of Eq. ( p^| ) is also included in Fig. || (b), which 
as expected becomes also indistinguishable from the rest 
of the curves. From the comparison of Eq. (|l2|) and Eq. 
(O), the scaling property can be written as, 



r]\ip tp (r]x,r]t)\ 2 = \i/> tp {x , t)\* 



(14) 



In what follows we shall apply Eqs. (|12| ) and ([bj]) to 
analyze the speed of the propagating pulse. In partic- 
ular, we contrast the propagation of the peak value of 
the probability density obtained by deriving it with re- 
spect to time, to the peak value obtained by deriving 
with respect to position. As we shall show below, the 
corresponding speeds differ by a factor of 3 1 / 2 . 

The maximum value of |^>t p | 2 at any given point can be 
obtained straightforwardly by deriving Eq. (|lj). How- 
ever, the following results show that the calculation of 
a maximum of \tptp\ 2 critically depends on whether the 
spatial or time evolution of the probability density is ana- 
lyzed. From (d/dt)\ip tp \ 2 = 0, we obtain the critical time 
which maximizes IV'tpl 2 at x/, the result is, 



V3' 



(15) 



where r = (xt /v qo ) . This result is in exact agreement 
with the scale found recently by Muga and Biittiker || 
for the time of arrival. 

It is clear that the pair (xf,t m ) gives the maximum 
value in a |^tp| 2 vs t plot. However, it is easy to show 
that the value xt does not correspond to the position of 
the peak of the traveling pulse. The latter, which is the 
maximum of a iV'tpl 2 vs x plot, is in fact located at the left 
of Xf (see below). From (d/ dx)\ip t p\ 2 = 0, at a fixed time 
t = tf, we obtain the maximum of the transient pulse as 
it propagates along the x coordinate, it becomes, 



= v an t 



qo L f- 



(16) 



Notice that for the particular case tf = t m , the maxi- 
mum of the pulse x m is retarded with respect to Xf i.e. 
x m = 3 _1 ^ 2 a;/. In other words, the maximum intensity 
of \ipt P \ 2 is attained at Xf, before the arrival of the pulse's 
peak takes place. Although the above statement seems to 
be in contradiction, we shall show below that these two 
types of maximum values have different meanings, and 
consequently there is no contradiction at all. The values 
of the probability density at these two positions are illus- 
trated on the \iptp\ 2 vs x plot (solid line) depicted in Fig. 
0; they are compared at the same time, t m = 30 fs. The 
values of \iptp\ 2 at x m and Xf are indicated on the curve 
by a full dot and a hollow circle, respectively. Equation 
(U|) shows that the maximum of the pulse (full dot) trav- 
els at exactly the semiclassical group velocity v qo , while 
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the other (hollow circle), according to Eq. (|l5|), travels 
faster than the former, at the speed 3 1 ^ 2 v qo . In fact, the 
observed spreading of the pulse is due to the fact that 
all the points of the IV'tpl 2 vs x curve are traveling with 
different velocities. 

It is clear from Fig. fj] that the height of the hollow 
circle hh c (t m ) is smaller than the corresponding height 
of the full dot, hfd(t m )- However, the height of the pulse 
diminishes at such a rate, that when the full dot reaches 
Xf at a later time t' = 3 1 / 2 t m , its height hfd(t') is smaller 
than hhcitm)- This is illustrated in Fig. p] where we plot 
|^tp| as a function of the position x, at a time t = t' 
(dashed line). The above effect can be easily shown by 
using the scaling property Eq. (|l4|). But first, lets us 
obtain the values of the heights hh c and h fd measured at 
t = t m , namely 

M*m) = \Tpt P (Xf,t m )\ 2 = ~ ( ^ — J ; (17) 

4 \2ir x m q J 

and 

hfd(t m ) = \i>t P (xm,t m )\ 2 = ( — — ] . (18) 

\Zirx m q Q J 

If we now feed r/ = 3 1 / 2 in Eq. (|lj), using the fact 
that t' — 3 1 ^ 2 t TO and xj = 3 1 ' 2 i m , we obtain that 
3 1/,2 /i/d(i') = hfd(t m ). This allows us to write the height 
ratio as, 

hhc{t m ) _ 3V3 , , 

h fd (t') 4 >L [U) 

Thus, we have in general demonstrated that the maxi- 
mum value of IV'tpl 2 ) measured at a fixed position Xf, is 
always reached there before the arrival of the maximum 
of the transient forerunner at that particular position. In 
view of this result, the role of the different time scales is 
now clear. While the traversal time r does in fact corre- 
spond to the time of arrival of the maximum of the tran- 
sient forerunner at Xf, the time scale 3 -1 / 2 r corresponds 
to the attainment of the maximum value of \iptp(xf, t)\ 2 
at Xf, which as demonstrated above, is always associated 
to values of the probability density located at the front 
tail of the forerunner. 

The above discussion shows that the calculation of the 
time of arrival may critically depend on whether the spa- 
tial or the time evolution of the probability density is 
analyzed. Although the interpretation of the results is 
different in each case, we believe that both kind of anal- 
ysis are important and complementary. It is important 
to point out that a time-frequency analysis of the fore- 
runner, performed in Ref . ^| , has revealed that it is com- 
posed predominantly by over-the-barrier (non-tunneling) 
frequencies. 

III. CONCLUSIONS 

The issue of spatial and time dependence of quan- 
tum waves in a step potential barrier has been analyzed 



4.0x10 




Distance (nm) 



FIG. 7: Spatial evolution of the transient pulse calculated 
from IV'tpi 2 for t m = 30.0 fs (solid line). Its values at the 
positions x m — 48.639 nm and xj = 84.246 nm are indicated 
respectively by a full dot and a hollow circle on the solid line. 
At a later time, t' = 3 1//2 i m , the maximum of the forerunner 
reaches the position Xf, as indicated by the full dot on the 
dashed line. Notice the reduction of its height 

by means of an exact solution to the time-dependent 
Schrodinger equation, for a point source initial condition. 
The main results of our study are the following: (1) For 
the case loq < V in the opaque barrier regime (xqo 3> 1), 
we found that the peak of the forerunner travels at ex- 
actly the group velocity v qo = (hq^/m), which implies 
that the time of arrival of the maximum of the peak at 
a fixed position in space Xf, is given by the traversal 
time t = (xf I v qo ) . This is in contrast to the time scale 
3 _1 / 2 t described in Ref. || . (2) We derived from the for- 
mal solution a closed analytical formula that accurately 
describes the dynamics of the forerunner, and found that 
it obeys an interesting and useful scaling property. (3) 
Although 3~ 1//2 r actually corresponds to the maximum 
of the forerunner when it is analyzed as a function of 
time (and fixed x), we demonstrated by means of the 
scaling property that this time scale always correspond 
to the passage of the front tail (not the maximum) of the 
forerunner across the fixed position Xf. 

As a final remark, it is important to point out that 
the forerunner is mainly composed by over-the-barrier 
(non-tunneling) frequencies, as recently shown by Muga 
and Biittiker ||. Thus, it is a surprising fact that the 
traversal time r, usually regarded as one of the possible 
tunneling times, is here associated to a non-tunneling 
process. 
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APPENDIX A: POINT SOURCE PROBLEM 

In this appendix we shall present a different derivation 
of the solutions to the point source problem 0,0], by 
means of the Laplace Transform method. Let us con- 
sider an exact analytical solution to the time-dependent 
Schrodinger equation for a step potential barrier V(x) — 
Q(x)V , 



2m dx' 



+ V(x)-ih 



i/j(x,t) = 0, 



(Al) 



subject to the source boundary condition given by Eq. 
(U). We must also assume that inside the barrier (x > 
0) the initial state satisfies the boundary condition 
ipo (x; * = 0) = 0, for x > 0. 

To obtain the solution for x > Oat t > 0, we begin by 
Laplace transforming equation ( |Al[ ) using the standard 
definition 



ijj(x;s)= / ip(x,t)e st dt, 



(A2) 



s = 7'. Aft er a simple partial fraction decomposition of 
in Eq. (A6) we obtain, 



2m 



7'+ioo ~y'-\-ioo 

.4.)*+ / •-<•>* 

7'— zoo 7' — 200 



> (A7) 



where the integrands $±(s) are defined as, 

e lpx e st 



$±(s) 



2 p(p±k )' 



(A8) 



In order to evaluate the integral (A7) we perform the 
change of v ariable, s' — s+iV, which allows us to identify 
in Eq. ( A7) the integral representation of the Moshinsky 
M -functions [@, 1 



M(x,±k ,t) 



7+ioo 

1 f i/3 



2m J 2 V^!?(vW±M 

7 — 200 



da', 



and obtain the solution for the case of u>q > V, given by 
Eq. (0). The solution for the case of luq < V, given by 
Eq. (H), is readily obtained along the same lines as Eq. 



with the initial condition given by Eq. (|l|). Let us first 
consider the case propagating waves (uq > V) which cor- 
responds to a dispersion relation fco = [/3(wo — ^)] 1 ^ 2 , 
where we have defined f3 = (2m/ h) and V — (Vo/h). 
The Laplace transformed solution reads, 



x>0, 



(A3) 



where p 2 = (3(is — V). The corresponding Laplace trans- 
form of Eq. (|l|) yields 



V (0; S ) 



t > 0. 



IUJ 



(A4) 



By combining Eqs. (A3) and (A4) evaluated at x — we 
can determine the value of the constant ci , which yields 
the Laplace Transformed solution for x > 0, 



ip(x;p) 



(p + ko)(p- h,)' 



x>0. 



(A5) 



The time dependent solution for x > is readily ob- 
taine d b y performing the inverse Laplace transform of 
Eq. (A5), using the Bromwich integral formula, 



tp(x,t) 



1 

2nd 



7 +zoo 



ip(x; s)e st ds, 



(A6) 



7'— zoo 



where the integration path is taken along a straight line 
s = j' parallel to the imaginary axis in the complex s- 
plane. The real parameter 7' can be chosen arbitrarily 
as long as all singularities remain to the left-hand side of 



1. Asymptotic behavior of the solutions 

In this subsection we consider the long time behavior 
of the solutions given by Eqs. (||) and (|[). This can 
be analyzed using the series expansion in powers of the 
argument y q (Eq. (§)) of the M-functions (Eq. (@)) [§. 
For large values of the argument y q i.e. \y q \ 1, the 
M -functions have the following series representation (see 
Ref. i): 



M(y q ) 



1 



1 /2htt 



1 



1 



n 1 / 2 y q n 1/2 yl 



(AlOa) 



provided that the phase, <fi q , of y q = \y q \ exp(0 9 ) lies in 

the interval — tt/2 < <fi q < ir/2. 

For the case tt/2 < <fi q < 3ir/2, instead of Eq. ( AlOa ) 
we have, 



M{y q )**\j mx 2 ' 2ht [2eyl + ^- 



n 1 / 2 y3 



(Alia) 



In general, the long time regime (t — > 00) corresponds 
to large values of the argument y q . For the particular 
values of q — ±iqo, we have that y q — > — e~ 47r / 4 (t//3) 1 / 2 g. 
Note that this expression depends on the choice of qo, 
and by inspection, we can see that for the cases q = iqo 
and q = —iqo the arguments are 4>i qo = — 57r/4 and 
A iq — ~ 7r /4, respectively. Therefore according to Eqs. 



(AlOa) and (Alia), one obtains that in the limit t — > 00, 
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M(y_ iqo ) -> and -> e -9oz e *(V-^ )t. Substitut- 

ing this values in Eq. (|j) yields the asymptotic solution 
given by Eq. (|^). By performing a similar analysis for 



the case of wave evolution for luq > V given by Eq. (||), 
one finds that the asymptotic behavior as t — > oo is given 
by Eq. (@). 
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